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Abstract. We discuss the statistical mechanics of rotating self-gravitating systems by allowing properly for the con- 
servation of angular momentum. We study analytically the case of slowly rotating isothermal spheres by expanding 
the solutions of the Boltzmann-Poisson equation in a series of Legendre polynomials, adapting the procedure in- 
troduced by Chandrasekhar (1933) for distorted polytropes. We show how the classical spiral of Lynden-Bell 
& Wood (1967) in the temperature-energy plane is deformed by rotation. We find that gravitational instability 
occurs sooner in the microcanonical ensemble and later in the canonical ensemble. According to standard turning 
point arguments, the onset of the collapse coincides with the minimum energy or minimum temperature state in 
the series of equilibria. Interestingly, it happens to be close to the point of maximum flattening. We determine 
analytically the generalization of the singular isothermal solution to the case of a slowly rotating configuration. 
We also consider slowly rotating configurations of the self-gravitating Fermi gas at non zero temperature. 
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1. Introduction 

Recently, the statistical mechanics of self-gravitating sys- 
tems has attracted considerable attention (Chavanis, 
Sommeria & Robert 1996, Chavanis & Sommeria 1998, 
de Vega, Sanchez & Combes 1998, Youngkins & Miller 

2000, Follana & Laliena 2000, Semelin, Sanchez & de Vega 

2001, Cerruti-Sola, Cipriani & Pettini 2001, Ispolatov & 
Cohen 2001, Chavanis 2002a, Taruya & Sakagami 2002, de 
Vega & Sanchez 2002, Huber & Pfenniger 2002...). This 
topic was introduced in the 1960s by Antonov (1962) and 
Lynden-Bell & Wood (1968) and further developed by 
Hertel & Thirring (1971), Horwitz & Katz (1978), Katz 
(1978) and Padmanabhan (1989) among others (see a 
complete list of references in the review of Padmanabhan, 
1990). These authors pointed out the particularity of 
self-gravitating systems to possess negative specific heats. 
They showed that this strange property is responsible for 
the inequivalence of statistical ensembles (microcanoni- 
cal/canonical) and the occurence of giant phase transi- 
tions associated with gravitational collapse. For a long 
time, these topics were only discussed in the astrophys- 
ical literature and were considered as a curiosity (not to 
say a fallacy) from statistical mechanicians. The situation 
is changing lately as these properties are re-discovered for 
other physical systems with long-range interactions which 
can be studied in the laboratory (see, e.g.. Gross 2001). 
For that reason, maybe, the statistical mechanics of self- 
gravitating systems comes back to fashion with new per- 
spectives. 



From our point of view, the statistical mechanics of 
self-gravitating systems is far from being completely un- 
derstood and rests on simplifying idealizations. The first 
idealization is to enclose the system within a box so as 
to prevent evaporation. It is only under this condition (or 
by introducing more realistic truncated models) that a 
rigorous statistical mechanics of self-gravitating systems 
can be carried out. Second, for most astrophysical sys- 
tems, the relaxation time by two-body encounters is much 
larger than the age of the universe so that a more subtle, 
colhsionless, relaxation must be advocated to explain the 
structure of galaxies. This is the concept of violent re- 
laxation formalized by Lynden-Bell in 1967. Then, it is 
implicitly assumed that the relaxation towards statistical 
equilibrium proceeds to completion, which is not necessar- 
ily the case in reality. Indeed, it is possible that the relax- 
ation stops before the maximum entropy state is attained 
(see Lynden-Bell 1967, Tremaine, Henon & Lynden-Bell 
1987, Chavanis et al. 1996). This problem of incomplete 
relaxation must be approached with extensive numerical 
simulations. Of course, this program was started long ago 
(e.g., van Albada 1982) but only recently are iV-body sim- 
ulations carefully compared with the predictions of the 
statistical mechanics approach (Cerriti-Sola et al. 2001, 
Huber & Pfeninger 2002) with variable success. 

The statistical equilibrium of a non-rotating classi- 
cal gas enclosed within a box was first investigated by 
Antonov (1962). He worked in the microcanonical en- 
semble and found that thermodynamical equilibrium ex- 
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ists only above a critical energy Emin = — 0.335GM^/i?. 
Below that energy, the system is expected to collapse 
and overheat; this is the so-called "gravothcrmal catastro- 
phe" (Lynden-Bell & Wood 1968). An isothermal collapse 
also occurs below a critical temperature T„j„ = 
in the canonical ensemble. This thermodynamical insta- 
bility is closely related to the dynamical Jeans instabil- 
ity (Semelin et al. 2001, Chavanis 2002a). In a recent 
scries of papers, we considered some extensions of the 
Antonov problem to the case of self-gravitating fermions 
(Chavanis & Sommeria 1998, Chavanis 2002c), special and 
general relativity (Chavanis 2002b) and confined poly- 
tropes (Chavanis 2002d). These studies complete previ- 
ous works on the subject. We also introduced a simple 
dynamical model of Brownian particles in gravitational 
interaction (Chavanis, Rosier & Sire 2002). By introduc- 
ing by hands a friction and a noise, we force the system to 
increase entropy continuously, thereby avoiding the prob- 
lem of incomplete relaxation. This model can be used to 
test precisely the ideas of equilibrium statistical mechan- 
ics (inequivalence of ensembles, phase transitions, gravita- 
tional instabilities...) and is sufficiently simple to allow for 
a thorough analjdiical investigation of the collapse regime 
when an equilibrium state does not exist. 

In all these studies, the system is assumed to be non- 
rotating so that the conservation of angular momentum 
is trivially satisfied. The object of the present paper is to 
extend the statistical mechanical approach to the case of 
rotating self-gravitating systems. This problem has been 
considered previously by Lagoute & Longaretti (1996), 
Lahena (1999), Lynden-Beh (2000) and Fliegans & Gross 
(2002) with different types of models. Clearly, the most 
interesting situation is the case of rapidly rotating sys- 
tems since a wide variety of structures can emerge as 
maximum entropy states (Votyakov et al. 2002). However, 
we shall restrict ourselves in the present paper to the 
case of slowly rotating systems for which the problem 
can be tackled analytically. We shall adapt to the case 
of isothca'mal spheres the classical proccdiire developed 
by Cliandrasckhar (1933) for distorted polytropes, i.e. we 
shall expand the solutions of the Boltzmann-Poisson equa- 
tion in terms of Lcgendre polynomials. A similar proce- 
dure was performed by Lagoute & Longaretti (1996) for 
rotating globular clusters subject to tidal forces and de- 
scribed by an extended Michic-King model. We believe 
that it is useful to consider the case of a gas enclosed 
within a box so as to make a clear connexion with the 
Antonov model when the rotation is set to zero. In par- 
ticular, we shall derive the expression of the thermody- 
namical parameters for slowly rotating isothermal spheres 
and show how the classical spiral of Lynden-Bell & Wood 
(1968) in the E — T plane is modified by rotation. We 
shall show that rotation avances the onset of gravothcrmal 
catastrophe in the microcanonical ensemble and delays the 
isothermal collapse in the canonical ensemble. Using the 
turning point criterion of Katz (1978), we argue that the 
series of equilibrium becomes unstable at the point of min- 
imum energy (in the microcanonical ensemble) or at the 



point of minimum temperature (in the canonical ensem- 
ble). Interestingly, these instabilities happen to be close, 
in each ensemble, to the point of maximum flattening. We 
establish the generalization of the singular isothermal so- 
lution to the case of a slowly rotating configuration. We 
also consider the case of slowly rotating self-gravitating 
fermions. This system exhibits phase transitions between 
"gaseous" states with an almost uniform distribution of 
matter and "condensed" states with a core-halo struc- 
ture. By cooling below a critical temperature, an almost 
nonrotating gaseous medium can collapse into a rotating 
"fermion ball" containing a large fraction of mass and an- 
gular momentum. 

2. Statistical mechanics of rotating 
self-gravitating systems 

2.1. The mean- field approach 

Consider a system of N particles, each of mass m, 
interacting via Newtonian gravity. We allow the sys- 
tem to have a non vanishing angular momentum. Let 
/(r,v,t) denote the distribution function of the system, 
i.e. /(r, V, f)(i^rd^v gives the mass of particles whose po- 
sition and velocity are in the cell (r, v; r -|- d^r, v -|- d^v) at 
time t. The integral of / over the velocity determines the 
spatial density 



/ 



(1) 



The mass and angular momentum of the configuration are 
given by 



M 



pdh, 



(2) 



(3) 



On the other hand, in the mean-field approximation, the 
energy can be expressed as 



L = / / r X V d^rd^v. 



E=^j fv^d^rd^v + p<^dh = K + 



(4) 



where K is the kinetic energy and W the potential energy. 
The gravitational potential $ is related to the star density 
by the Newton- Poisson equation 



A$ = AwGp. 



(5) 



The equilibrium configuration of the system is deter- 
mined by max;imizing the Boltzmann entropy 



-k 



[liUd^ 

J m m 



rd^v, 



(6) 



while conserving mass, angular momentum and energy. 
We stress that for systems interacting via a long-range po- 
tential, like gravity, the mean-field approximation is exact 
so that this procedure is entirely rigorous and provides a 
simple approach to the problem. 
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Introducing Lagrange multipliers a, (3 and — /30 for 
each constraint, we find that the critical points of entropy 
are given by 



/ = ^'g-/3(^ + *-S^-(rxv))^ 



where 
m 



(8) 



We can rewrite the foregoing expression for / in the more 
suggestive form 



-/3(<i>-i(Oxr)") 



(9) 



We note that the most probable distribution of a rotating 
self-gravitating system is a Maxwell-Boltzmann law with 
a dispersion the same at every point (isothermal distribu- 
tion). Moreover, the most probable form of rotation is a 
rigid rotation Q,. By integrating Eq. (||) over the velocity, 
we find that the density is given by the Boltzmann factor 



(10) 



The quantity in parenthesis is the "effective" potential in 
the rotating frame 



$e//-$-i(f2xr)2. 



(11) 



accounting for inertial forces. The equilibrium configura- 
tion is obtained by solving the Boltzmann-Poisson system 



(12) 



and relating the Lagrange multipliers A, (3 and SI to 
the constraints M, E and L. The configuration is ther- 
modynamically stable if the second order variations of 
entropy are negative for any perturbation that satisfies 
the constraints on mass, energy and angular momen- 
tum. This corresponds to the microcanonical description. 
Alternatively, we could describe the system in the canoni- 
cal ensemble assuming that (3 and are given, instead of 
E and L. In that case, the equilibrium configuration is ob- 
tained by maximizing the free energy J = S — j3E — /3r2 • L 
at fixed M, /3 and fJ. Of course, the microcanonical and 
canonical ensembles yield the same critical points, i.e. the 
critical points of entropy at fixed mass, energy and an- 
gular momentum and the critical points of free energy 
at fixed mass, temperature and angular velocity coincide. 
Only the onset of instability, regarding the second order 
variations of S or J, will differ from an ensemble to the 
other. In the non-rotating case, the thermodynamical sta- 
bility analysis was performed by Padmanabhan (1989) in 
the microcanonical ensemble and by Chavanis (2002a) in 
the canonical ensemble, by solving an eigenvalue equa- 
tion. The onset of instability can also be determined by 
the turning point criterion of Katz (1978) who extended 
the theory of Poincare on linear series of equilibrium. 



2.2. The rotating isothermal sphere 

To determine the structure of rotating isothermal spheres, 
we first introduce the function 5* = (3{^eff — ^o), where 
$0 is the gravitational potential at r = 0. Then, the den- 
sity field can be written 



P = Poe , 



(13) 



where pQ is the central density. Introducing the notations 
^ = {i-KCPpoY^^r, V = Q^/2'kGpo and using spherical 
coordinates r, 9, (j), we can rewrite the Boltzmann-Poisson 
equation (113) in the form 



(14) 



where p, — cos 9 and we have neglected (j) by assuming 
that the configuration is symmetric with respect to the 
axis ft = riz. The true gravitational potential <f> is related 
to ^' by the relation 



/3$ = ^ + ^vC'^{l^ p^)+ p<i>o. 



(15) 



Eq. ( p^ is the fundamental equation of the problem. For 
V ^ (no rotation), we recover the Emden equation 



_ ,2^ 



-V 



(16) 



with ■(/) = '(/)' = Oat^ = (Chandrasekhar 1942). 

So far, we have made no approximation regarding the 
value of the angular velocity. We shall now consider the 
case of slowly rotating structures and let w — > 0. Assuming 
the following form for our solution 



■ii ^ ifj + vQ + v^Q2 



(17) 



and keeping terms only to first order in u, we find that O 
satisfies the equation 

(18) 

Now, following a procedure that dates back to Milne 
(1923) and Chandrasekhar (1933), we shall assume for 
the following form 



e = 0o(O 



+ 00 



^A,4>,{£,)P,{p), 



(19) 



where Pj{p) are the Legendre polynomials satisfying the 
differential equation 



d 



dp 



dp 



— (l-M^)^ +,7(.? + l)P,=0. 



(20) 



Substituting for O from Eq. (^9|) in Eq. (|T^) and equating 
coefficients of Pj, we get 



1 d / 2d(f>o 
ed^V d^ 



(21) 
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1 d ^^2'^'t']\ I i(j + 1) ^'^'^0 '^^^ ^^'-'^^y ^"^S^l^r isothermal sphere 

^ d^K J \ i J Pqp ^ ^ +00, the solution of the Emden equation (16) 

{j = 1, 2, ...) (22) behaves like (Chandrasekhar 1942) 



with <?!)j = 0^. = at ^ = 0. (32) 

So far, the Aj are arbitrary. They will be determined 

by requiring that the gravitational potential and its radial Substituting this asymptotic behaviour in Eq. (|^) and 

derivative are continuous across the sphere at r = i?. Now, making the change of variables ^ = e*, we obtain 

outside the sphere the potential is given by the Laplace ^ 

equation + + 2<^o = -e^*. (33) 



A$e.t = 0. (23) 



dt^ dt 

The solution of this equation is readily found. After re- 



rriT. u • n i 1.1 1 i- f ii • i J turning to original variables, we get for large f 's: 

The physically acceptable solution of this equation tend- as j & as 

ing to zero at infinity can be written as 14/ ^/7 \ 

P^ext — ^ ~^ '^''^^ (^^) Following the same procedure for the function (^2 , we have 

I — k d(j)2 

On, the other hand, according to Eq. (|15D, we have inside — j-^ — I — — 4(p2 = 0, (35) 



the spherical box 
1 



di2 dt 
leading to 



/3$„* = * + |<2(i_p^(^))_^^$^^ (25) ^, 



cosh(^ln^ + 5'j. (36) 
where ^'(^,/i) is given by Eqs. (p7|) (pj|) and use has been 

made of the identity 1 - /i^ = |(1 - P2(m))- Let us denote Keeping only the leading term in Eqs. (||)(|6|) in the 

by a = {^TrGPpay^^R the value of ^ at r = i?. Comparing limit ^ +oo, we find that 
the inner and the external potentials at ^ = a, and also 

their derivative, we find that Aj ~ Bj ^ if j ^ 2 and, ^ — -C^; 02 K^'' , 7 = — -. (37) 

for j = 2, 8 ' '22 

^ Substituting these results in Eq. ( p9| ) and returning to 

= A2(t)2{a) - -t'' 
a-' 6 



^ — ^202 (a) o? , (26) dimensional variables, we finally obtain 



1 r 20 /i?^'"^ 



-^=^20^(«)-ia. (27) ^=2^ + 8;?Gr-3(3T^bj ^2(^^)1' (^8) 

which is the expression of the slowly rotating singular 
Solving for A2, we get • 11, 

^ ' ° isothermal sphere. 

A2 - ^ (28) 

6 302(a) + "02 (o^) 3. The thermodynamical parameters 

Therefore, the solution of the Boltzmann-Poisson equation 3.1. The mass 
(fl4) to first order in v is given by 

In spherical coordinates, the mass is given by 

-^ = -^(^)+-{^°(^)n 302(.)+U(a) ^-^^)^-(^)}' M.2. r Tp^^d.^,. (39) 

(29) ^" 

Introducing the dimensionless variables previously de- 
with fined, we obtain, using Eq. (|13[), 



i? 2Q! 



-1 Jo 



Substituting for 4" from Eq. (^) in Eq. ( [io| ) and recalling 
^ A ^^2^^ ^ _ g-^j (3-,^^ that j;^^ Pj{lj)d^i = 0, we obtain to first order in v. 

with 0, = = at C = 0. f ^""(1 - -'Ao(0)e^de (41) 



Now, 
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3.3. The energy 



(42) 



a 



a 



= -y -a^(/>o(a). 



(43) 



Therefore, the normalized inverse temperature rj is ex- 
pressed in terms of a and v by the relation 



77 = aip'{a) + V 



(44) 



Expressing the central density in terms of a, we find that 



277 



with 



GM ) 



(45) 



Therefore, the foregoing expression for 77 can be rewritten 
to order cj^: 



6 a 



(46) 



The validity of our pcrturbativc approach requires that 



2w" 



^ + -'/'o(a) 



< 1. 



(47) 



3.2. The angular momentum 



Using the MaxwcU-Boltzmann distribution (^), we can 
rewrite the angular momentum (|^) in the form 



L = p r X u (fir, with u = x r. 



(48) 



After straightforward manipulations, we find that the an- 
gular momentum is related to the angular velocity by 



(49) 

)cPr is the moment of inertia 
tensor. If the density depends only on r and 9, we have 
L = /r2 with 



where Itj = J p{r^S: 



-1 pR 

/ = 27r / / pr'^il - p^)drdp. 
-1 Jo 



(50) 



For our purpose, it is sufficient to determine the angular 
momentum to first order in 17. Therefore, we need just the 
expression of the moment of inertia / for a non-rotating 
isothermal sphere. Introducing the dimensionless variables 
previously defined, and using the expression ( ^ ) for 77, we 
get 



Mi?2 3Q4^'(a) 7o 



Then, 



A = T{a)uj, with A 



L 

VGAPR 



(51) 



(52) 



Quite generally, the potential energy of a self-gravitating 
system can be written in the form (see, e.g., Binney & 
Tremaine 1987) 



W 



p r • V$ cfr. 



(53) 



It can be readily verified that the condition f = + 
<I>e//), where w = v — u and $e// is the effective po- 
tential (^1]), is equivalent to the condition of hydrostatic 
equilibrium in the rotating frame 



Vp = -/9V$ - pflx {nx r), 
where 



(54) 



(55) 



is the local pressure. For an isothermal gas p = J^pT. 
Inserting Eq. (^) in Eq. (p3|), we obtain after straightfor- 
ward calculations 



T^-t-3 / pd-^r + L-n - *pr-dS = 0, 



(56) 



where dS is a surface element normal to the spherical box. 
On the other hand, we can write the kinetic energy in the 
form 



K = Krot+Kth, 

where 

Krot = ^ y pu^d^r = ^I.-n, 
is the rotational energy and 
Kth^l I fw^dhd^w = 11 pd^r, 



(57) 
(58) 

(59) 



is the thermal energy. For an isothermal gas Kth — ^NkT. 
Therefore, Eq. (M) becomes 



W + 2K= (hpr-dS, 



(60) 



which is the Virial theorem for a self-gravitating gas en- 
closed within a box. For the total energy E = K + W, we 
have 



E = -K 



pr-dS. 



(61) 



This expression is valid even if the system is not axisym- 
metric. If, now, the density depends only on r and 9, and 
if the gas is isothermal, it is possible to combine the for- 
mulae derived in Sect. 2.2 to obtain the expression 



A; 



ER 
GAP 



3 

2^ 



1 



-Xlo 



-'^("'(l-«0o(a)), (62) 



where 77 and v can be expressed in terms of a and u by Eqs. 
(p|)(|45l). We also recall that our theory is valid to order 



For Lo = 0, Eqs. (Wq)(p3) reduce to the equations of 



state obtained by Lynden-Bell 
rotating isothermal spheres. 



Wood (1968) for non- 
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3.4. The entropy 

Using Eqs. @(|lO|) and ([l3|), the distribution function of 
a rotating isothermal gas can be written 



3/2 



(63) 



Substituting this expression in Eq. (||) and expressing po 
in terms of a, we get 



mS 



1 



-Mln(3 -2Mlna 



(64) 



Throughout this paper, we shall not write the constant 
terms (depending on the fixed parameters M and R) 
which enter in the expression of the entropy. Therefore, 
a term (3Kth = has been ignored in Eq. (64). Using 



the definition of ^P, the last integral can be rewritten 



Using W = E — Krot — Kth, we obtain 

^=-iA/ln/3-^/3L.ll 

-2M In a + 2(3E - (3M^a. 



(65) 



(66) 



Introducing the dimensionless parameters previously de- 
fined, we get 

^ = -i Inr; - ^77Ac^ - 21na - 277A - /3$o- (67) 

We now need to determine the central potential $0. The 
condition that $ and d^/d£, are continuous at ^ = a im- 
plies for j = 0: 



13 1 
a 6 



13 1 

— = ip'{a) + v(j)'^{a) + -va. 



(68) 



(69) 



Comparing Eq. ( |69[ ) with Eq. (|44D , we see that —Bq/u = 
77. Inserting this result in Eq. (|68|), we obtain 



— + (/)o(a) 
6 



(70) 



With this new relation, the entropy (|6^) becomes 
S 1 



■ In 77 — 2 In a -|- ^{a) + rj — 2A77 



--rjXuj + V 



— -I- 00 (a) 
6 



(71) 



3.5. The flattening function 



The value of the potential at ^ for a non-rotating configu- 
ration is ipiO- 9, rotating configuration, the equation 



of the surface with the same value of the potential is given, 
to first order in u, by 



V^'(0 



MO 



6 302 (a) + a^jlo^) 



02(O^'2(m) 



(72) 



If a denotes the largest radius of the isodensity surface (at 
the equator p = 0) and b the smallest radius (at the pole 
p = 1), and if we define the flattening by/ = l — 6/a (see, 
e.g., Lagoute & Longaretti 1996), we get 

^ n. „,2 



/«) 



4CV'(0 302(a) + a02 (a) 



MO- 



(73) 



Considering the limit ^ — s- +00 and returning to dimen- 
sional variables, we find that the fiattening function of the 
rotating singular isothermal sphere is 



fir) 



5 pn^R"^ 



R 



(74) 
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a result which can also be derived directly from Eq. (p8[). 

Coming back to Eq. ([73|), we find that the flattening 
behaves with the distance (for a given value of a) as 



(75) 



On the other hand, the flattening at the edge of the con- 
figuration (^ = a) is given by 



/(«) 



02 (a) 



(76) 



2 302 (ck) + Q;02(a) ' 

where we have used Eqs. (^5|)(^6|) to eliminate the variable 
V in profit of uj. In the canonical ensemble (fixed w), this 
function depends on a as 

02(a) 



(77) 



302(a) -I- a02(a) 

In the microcanonical ensemble (fixed A), we must express 
the angular velocity in terms of the angular momentum 
using the relation (^2|). In that case, the flattening at the 
edge of the configuration depends on a as 

FMCEia) = ^^^FcE{a). (78) 

We shall come back to these results in the following sec- 
tion. 

4. Equilibrium phase diagram 

4.1. Microcanonical ensemble 

The microcanonical ensemble corresponds to isolated sys- 
tems characterized by their energy A and their angular 
momentum A. In order to determine the equilibrium phase 
diagram 77(A) for different values of A, we need to solve 
Eqs. (16)(|30|)(|3l|) numerically. Expanding the functions 
'0, 00 and 02 in Taylor series for ^ ^ 0, we get 



^ 6^ 



120 



(79) 
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Fig. 1. Moment of inertia for a non-rotating isothermal 
sphere along the series of equilibria (parametrized by a). 



1 



1 

120^ 



(80) 



(81) 



so that V"(0) = 1/3, (t)'^{0) = -1/3 and </>^'(0) -2. The 
integration can be continued numerically by a standard 
Runge-Kutta routine. In Figs. |l|,||,|[ we plot the curves 
X(a), A(a) and riia) defined by Eqs. (|5l]) (H) (|^) . These 
curves exhibit damped oscillations and tend to the values 



= g, = 2 



— A^, As = A^. 

4 ' 4 32 



(82) 



as a — > +00. This asymptotic limit corresponds to the 
singular solution (|3^) . The iso-density contours of the ro- 
tating singular isothermal sphere are represented in Fig. 
I (for A = 0.17). 

In Fig. H, we have represented the curve ry(A). It has 
a classical spiral behaviour as noted by a number of au- 
thors in the non-rotating case. There is no equilibrium 
state (i.e., no critical point of entropy) above the value 
Ac(A). In that case, the system will collapse and overheat 
(gravothermal catastrophe) . It is also at this point that the 
critical points of entropy become unstable (saddle points) 
in the series of equilibria (Katz 1978). We see that rota- 
tion tends to favour the instability, i.e., the gravothermal 
catastrophe occurs sooner than in the non-rotating case. 

In Fig. 1^, we have represented the angular velocity of 
the system as a function of energy for different values of 
the angular momentum. We observe that the curve has a 
spiral behavior similar to the 77 — A diagram, but reversed. 
We might expect that the moment of inertia decreases 
as the system becomes more and more concentrated, re- 
sulting in an increase of angular velocity along the series 
of equilibria. This is true for moderate density contrasts 
(up to ~ 709), coinciding with the region of stability, but 



Fig. 2. Normalized energy of an isothermal sphere along 
the series of equilibria (parametrized by a) for different 
values of the angular momentum A = 0,0.1,0.2,0.3. 
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Fig. 3. Normalized inverse temperature of an isother- 
mal sphere along the series of equilibria (parametrized 
by a) for different values of the angular momentum A = 
0,0.1,0.2,0.3. 



not for larger density contrasts. Indeed, although the cen- 
tral density tends to diverge, the mass contained in the 
core is low and does not dominate the moment of inertia. 
Therefore, the moment of inertia and the angular veloc- 
ity have a non monotonous (in fact oscillatory) behaviour 
with a and saturate to finite values X, and ujg = X/Tg as 
a — > +00. 

In Fig. ^, we plot the flattening function F(^) defined 
by Eq. (|75|). As expected, the flattening is is a monotonous 
function of the distance. For f ^ 0, F ^ 3 and for ^ 
+00, F ~ In Fig. Ill, we plot the flattening at the 
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Fig. 5. Equilibrium phase diagram of isothermal spheres 
giving the inverse temperature r; as a function of minus the 
energy A for different values of the angular momentum A. 
The gravothermal instability at Ac occurs sonner (i.e. for 
larger energies) when the system is rotating. 



edge of the configuration Fmce{o) as a function of a in 
the microcanonical ensemble (see Eq. (^8|)). We observe 
that the curve displays damped oscillations towards the 
asymptotic value 3:^^(9/2)^. In particular the flattening 
(by unit of A^) is maximum for a — 32.6.... Interestingly, 
this value lies precisely in the range of values at which the 
gravothermal catastrophe sets in (compare with Fig. H). 

4.2. Canonical ensemble 

The canonical ensemble is characterized by the specifica- 
tion of the inverse temperature 77 and the angular velocity 



Fig. 7. Spatial dependance of the flattening function F{£). 



io. In Fig. ^, we have represented the curve A(?7) for dif- 
ferent values of w. There is no equilibrium state (i.e., no 
critical point of free energy) above the value ?7c(w). In that 
case, the system will undergo an isothermal collapse. It 
is also at this point that the solutions become unstable 
(saddle points of free energy) in the series of equilibria 
(Katz 1978). We see that rotation tends to delay the in- 
stability, i.e., the isothermal collapse occurs later than in 
the non-rotating case. 

In Fig. ^ we plot the flattening at the edge of the 
configuration Fq e (o^) as a function of a in the canonical 
ensemble. The curve displays damped oscillations towards 
the value The flattening (by unit of lo^) is maxi- 

mum for a = 5.4.... Interestingly, this value is close to 
the typical values at which the isothermal collapse sets 
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Fig. 8. Flattening at the edge of the configuration (per 
unit angular momentum squared) along the series of equi- 
libria (microcanonical ensemble). The flattening is maxi- 
mum for a — 32.6. This value lies typically in the region 
where the gravothcrmal catastrophe sets in. 
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Fig. 9. Equilibrium phase diagram of isothermal spheres 
giving minus the energy A as a function of the inverse 
temperature rj for different values of angular velocity. The 
isothermal collapse at ric{uj) occurs later (i.e. for smaller 
temperatures) when the system is rotating. 



Fig. 10. Flattening at the edge of the configuration (per 
unit of angular velocity squared) along the series of equi- 
libria (canonical ensemble). The flattening is maximum 
for a = 5.4. This value lies typically in the region where 
the isothermal collapse sets in. 




Fig. 11. Normalized inverse temperature of an isothermal 
sphere along the series of equilibria (parametrized by a) 
for different values of angular velocity a; = 0, 0.4, 0.8. The 
series of equilibria becomes unstable after the first maxi- 
mum. 



in (compare with Fig. |ri| giving the inverse temperature 
along the series of equilibria). 

5. Rotating self-gravitating fermions 

The previous results can be easily generalized to the case 
of self-gravitating fermions. This extension is relatively 
straightforward and we shall just give the main steps of 



the calculations. The thermodynamical parameters of a 
non-rotating Fermi gas at finite temperature have been 
calculated by Chavanis & Sommeria (1998) and we shall 
adopt a similar presentation. For a rotating configuration, 
the Fermi-Dirac distribution can be written 
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where 779 is the maximum allowable value of the distri- 
bution function and A > a strictly positive parame- 
ter (inverse of the fugacity) insuring that f < tjq (the 
other quantities have been defined previously). For quan- 
tum particles with spin s, rjQ = (2s -|- l)m'*/(27r?i)^. The 
distribution function ( p3| ) can be obtained by maximizing 
the Fermi-Dirac entropy 

5 = - ^l^\n^+(l-^)\n(l-^)}d^rd^^r,{84) 

J I '70 m V vaj V mj) 

at fixed mass, energy and angular momentum. This form 
of entropy also occurs in the context of violent relax- 
ation for coUisionless self-gravitating systems (Lynden- 
Bell 1967, Chavanis, Robert & Sommeria 1996, Chavanis 

6 Sommeria 1998). Defining * = Pi^eff - *o) and 
k = Ae''*" , the spatial density can be written 

p3/2 h/2{k), (85) 

where J1/2 is the Fermi integral 



In{t) 



-{-OQ 



-dx 



(86) 



of order n = 1/2. Substituting the relation (|85| ) in the 
Poisson equation (|^) and introducing the dimensionless 
parameters ^ = (47rG/3po)^^'^'' and 1; 



2irGpo 



we obtain 



^d_( ^d^\ ld_( ^d^\ _ h/2ike^) 

(87) 

We can also check that Eq. (|l^) keeps the same form. For 
u = 0, Eq. (|^ reduces to 



1 d 



ed^v d^ 



h/2{k) ' 



with ip = ip' = at = 0. For k +00 (classical limit), 
we can use the limiting form of the Fermi integral 



In{t) ~ |r(n + l) 



(t ~* +00) 



(89) 



and we recover the classical Emden equation (16). For 
k (completely degenerate limit), we have 



In{t) - ^ \ , (t^O) 



1 



(90) 



and the differential equation ( pq ) becomes equivalent to 
the Lane-Emden equation for a polytrope of index 3/2. 

We shall now consider the case of slowly rotatin g st ruc- 
tures and let 1; ^ 0. Repeating the steps of Sec. ^.2| , we 
find that Eq. (|29| ) remains valid with the new functions 
(pa and (f>2 defined by 



1 d 
^dl 



di 



ll_i/2{ke 
2 



h/2{k) 



1, 



(91) 



l/-i/2(fce'^) 



2 h,2{k) 



(92) 



with = at ^ = 0. In arriving at (|9|)-(|92|), we 

have used the identity 



-In-lit), {n>0), 



(93) 



for n — 1/2, which can be easily established from Eq. (^6|). 

We now determine the thermodynamical parameters 
of a slowly rotating self-gravitating Fermi gas. We can 
check that the inverse temperature is still given by Eq. 
( ^6| ) and that the relation (|5|) remains valid. On the other 
hand, eliminating the central density between Eq. (^) and 
a = {AirCPpof^R, we find that 



77 



^,h,2{kf 



where ^ rioV^Vh^Hi^MR^ , (94) 



is the "degeneracy parameter" (Chavanis & Sommeria 
1998). Using Eq. (||), we obtain 



/^'/i/2(fc) =aV(«)0 + 2c^ 



\ + -00 (a) 
3 a 



(95) 



The relation between the angular momentum and the an- 
gular velocity is still given by Eq. (^ ) with the moment 
of inertia 



3a V (a) h hMk) 
Using Eq. (|6l]), the energy can be put in the form 

a' r(h/2{ke^) 3 h/2{ke^) 



(96) 



A 



Jo 



/l/2(fc)5 2 /i/2(fc)^ 



MO >Cd^ 



1, 2aio r/3/2(Ae'^(")) 3 /i/2(fce'^("))_^ ^ ^ 
2 3 fJ.^ I h/2{kf 2 h/2{kf 



Finally, for the entropy we obtain (see Appendix ^) 



(97) 



VoS , , , , , 7 , 4 

— = In + 77 + V(a) - g - -ilXuj + ■ 



-7- + 'Po(a) 



2 

9^ 



/3/2(fce'/'(")) 3 /i/2(A:e'^(")) 



A/2(fc)= 



(98) 



For w = 0, Eqs. (||) (|9|)(£3) and (g8|) reduce to the equa- 
tions of state obtained by Chavanis & Sommeria (1998) for 
the non-rotating Fermi gas. For k — > -f 00 (non degenerate 
limit), we recover the equations derived in Sec. ^. 

The equilibrium phase diagram can be obtained in the 
following manner. For given fc, and uj, we can solve Eqs. 
(H) (|l|) and (H) until the value = a for which the 
relation (|5|) is satisfied. Then, Eqs. (^) and (|9^) deter- 
mine the temperature and the energy of the configuration. 
If the angular momentum is fixed instead of the angular 
velocity, we must use Eq. ( p^ with Eq. (^ to express 
Lo in terms of A. By varying the parameter k (for a fixed 
value of the degeneracy parameter /x), we can cover the 
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Fig. 12. Equilibrium phase diagram of self-gravitating 
fermions giving the inverse temperature 77 as a function 
of minus the energy A for different values of the angular 
momentum A and for a degeneracy parameter fj, = 10^. 



whole diagram in parameter space. A complete descrip- 
tion of this diagram has been given by Chavanis (2002c) 
in the non-rotating case. 

In Fig. ^ we represent the equilibrium phase diagram 
of self-gravitating fermions for a degeneracy parameter 
/i = 10^ and for different values of angular momentum. 
We observe that degeneracy has the effect of unwinding 
the spiral of Fig. |^. For sufficiently large values of the 
degeneracy parameter, there is still gravitational collapse 
at Ac accompanied by a rise of temperature, but this 
"gravothermal catastrophe" stops when the core of the 
system becomes degenerate. This leads to the formation 
of a "fermion ball" which contains a moderately large frac- 
tion of mass aM (at point D, we typically have a ~ 0.2). 
In the microcanonical ensemble, the decrease of potential 
energy in the core is compensated by an increase of tem- 
perature. Therefore, the mass (1 — a)M contained in the 
halo undergoes an expansion which, in our model, is ar- 
rested by the walls of the box. As a result, the density 
of the halo is almost uniform. Typical density profiles are 
given by Chavanis & Sommeria (1998) in the non rotating 
case. The expansion of the halo explains why the moment 
of inertia of the system increases despite the formation 
of a massive nucleus. Therefore, the angular velocity de- 
creases during the collapse contrary to what might be ex- 
pected. The angular velocity is represented as a function 
of energy in Fig. |l^. It has a complicated behaviour which 
corresponds to the unwindement of the spiral of Fig. ^. 

For smaller values of the degeneracy parameter, the 
gravitational phase transition is suppressed (Chavanis & 
Sommeria 1998, Chavanis 2002c) and the equilibrium 
phase diagram has the structure of Fig. 14, This diagram 
is similar to the one found by Fliegans 



Gross (2002) in 
their two-dimensional model of rotating self-gravitating 



-0.4 -0.1 
A=-ER/GM' 

Fig. 13. Angular velocity vs energy for self-gravitating 
fermions with a degeneracy parameter ij, = 10^ and an 
angular momentum A = 0.3. 

systems. Had they used a smaller value of their cut-off 
radius (which plays the role of the inverse of our degen- 
eracy parameter), they would have probably obtained a 
diagram similar to the one of Fig. However, when the 
cut-off radius (or degeneracy) is sufficiently large, the spi- 
ral unwinds and the 77- A curve is univalued like in Fig. 
|l|. For low energies, the equilibrium states have a core- 
halo structure with a partially degenerate nucleus and a 
dilute envelope. As energy decreases any further, the nu- 
cleus contains more and more mass and becomes smaller 
and smaller (a property of the R ~ M~^^^ law of degen- 
erate configurations). For rotating systems described in 
the microcanonical ensemble (fixed E and L), the result- 
ing decrease of moment of inertia is accompanied by an 
increase of angular velocity as shown in Fig. In this 
diagram, the uj — A spiral of Fig. ^ is completely unwound 
and the angular velocity increases monotonically with A. 

In Fig. ^ we have focused our attention to what hap- 
pens close to the minimum energy. For E — s- Emim T 
and the system has the same structure as a cold white 
dwarf star (see Appendix For our perturbative anal- 
ysis to be valid, we have taken a very small angular ve- 
locity A — 0.01. We shall describe the diagram of Fig. [l| 
in the {L,T) ensemble in which the angular momentum 
and the temperature are assumed given (this is a situa- 
tion intermediate between microcanonical and canonical 
ensembles). For T > Tc, the system is in a gaseous phase 
with a smoothly decreasing density profile. For T < Tc, the 
system undergoes an isothermal collapse that only stops 
when gravity is balanced by the degeneracy pressure. The 
result of this phase transition is a "fermion ball" which 
contains almost all the mass, unlike in the microcanonical 
ensemble at the point of gravothermal catastrophe. For 
rotating systems, this isothermal collapse is accompanied 
by a discontinuous rise of angular velocity at the critical 
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Fig. 14. Same as Fig. [ij for /i = 10'^. For high degen- 
eracy, the spiral is unwound and the gravitational phase 
transition is suppressed in the microcanonical ensemble. 
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Fig. 15. Same as Fig. y for ^ = 10^. 

temperature Tc (see Fig. |l^). Therefore, even if the initial 
rotation of the system is negligible in the gaseous phase, 
after collapse the "fermion ball" can have appreciable ro- 
tation as suggested in Fig. Its structure is then similar 
to a distorted polytrope of index n = 3/2 as computed by 
Chandrasekhar (1933) in the limit of slow rotation. 

6. Conclusion 

In this paper, we have considered the effect of a small ro- 
tation on the thermodynamical stability of self-gravitating 
systems. We have worked in a finite box in order to make 
a clear connexion with the Antonov problem for non- 
rotating systems and render the statistical mechanics of 
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Fig. 16. Equilibrium phase diagram for self-gravitating 
fermions with a degeneracy parameter fi = 10^ and an 
angular momentum A = 0.01. For T < Tc, i.e., 77 > rjc, 
the system undergoes an isothermal collapse leading to a 
rotating fermion ball containg a large fraction of mass and 
angular momentum. 



these objects rigorous. Physically, this idealization means 
that our isothermal system is surrounded by a medium 
which fixes its size. We have found that for rotating sys- 
tems, the well-known inequivalence of statistical ensem- 
bles for self-gravitating systems manifests itself in a strik- 
ing manner: the instability is advanced in the microcanoni- 
cal ensemble and delayed in the canonical one. In addition, 
we have found a connexion between the onset of instability 
and the configuration of maximum flattening in the series 
of equilibria. These results have be generalized to the case 
of self-gravitating fermions. 

On the other hand, the case of rapidly rotating isother- 
mal configurations is interesting because new, non trivial, 
structures can emerge as maximum entropy states (i.e., 
most probable states). The classification of such struc- 
tures is complicated because many bifurcations can oc- 
cur depending on the values of the control parameters 
(M, E, L). In particular, we must be careful to select only 
entropy maxima^ discarding the critical points of entropy 
which are only saddle points. This can be obtained either 
by solving the Boltzmann-Poisson equation (Votyakov et 
al. 2002) and checking the stability of the solutions or by 
using relaxation equations towards the maximum entropy 
state (Chavanis et al. 1996). These studies are important 
in order to obtain a classification of the most probable 
configurations of self-gravitating systems. The second ap- 
proach (relaxation methods) will be explored in a future 
paper (Rosier & Chavanis, in preparation). 
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Fig. 17. Angular velocity vs temperature plot for /i = 
10^ and A — 0.01. The gravitational collapse below 
is accompanied by a discontinuous jump of angular ve- 
locity due to the moment of inertia decrease. The re- 
sulting fcrmion ball rotates with angular velocity ~ 
0.1891Gi/2Mi/2i?-3/2 (see Appendix ||) corresponding 
to an increase ^ 270 of the initial angular velocity. The 
dimensionless rotation rate v ~ 3.98 10~^ is sufficiently 
small to justify the perturbative approach followed in that 
paper. 
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Appendix A: The entropy of a rotating 
self-gravitating Fermi gas 

The Fermi-Dirac entropy (|8^ ) can be written in the equiv- 
alent form 



/ 



/ 



ln( 1 - — 1 \Sxd^- 



Inserting the distribution function 



f = 



1 + fce^e^"^ 



(A.2) 



in Eq. (AT), we obtain 



-qoS = M In fc + / p^d^r + (3Kth 



-7/0 / ln( 1 + ■^e"*e-'^'f ]d^rd\. 



(A.3) 



The last integral can be reduced by an integration by parts 
to the expression {2/?>)f3Kth- Therefore, 



77o5 = MlnA: + / p^d^v + -(3Kth. 



(A.4) 



Tl^S = Aflnfc + 2f3W - PKrot - M/3*o + ^PKth- (A.5) 

Since W ^ E — K and Kth — K — Krot, the foregoing 
expression can be rewritten 

T]aS = Mlnfc + 2(3E - ^PK - Af/3$o - ^pK^ou (A.6) 
or, according to Eqs. ( p8| ) and (|6l]), 

r]oS = Mlnfc ^PE - i/3 (f prdS - Af/3$o - ^/3L • fl. 

o o J o 

(A.7) 

The pressure at the surfa ce of the sphere can be calculated 
with Eqs. (55) and (A^). After simplification, we get 



M = lnfc-I^A-/3<i>o-^?7Ao. 



h/2{k) 



5 2"" /i/2(fc)5 



Ma)}- (A. 



Using Eq. (|70D, which remains valid for self-gravitating 
fermions, and Eq. (^) we finally arrive at the form (|9^). 

Appendix B: The energy of a rotating polytrope 

In this Appendix, we derive a simple analytic formula for 
the potential energy of a slowly rotating polytrope of index 
n. The index rt = 3/2 describes a completely degenerate 
Fermi gas at zero temperature, which is a particular limit 
of the model studied in Sec. |^. 

For a self-gravitating system rotating with constant 
angular velocity il, the condition of hydrostatic equilib- 
rium in the rotating frame can be written 



(B.l) 



where ^e// is the effective potential defined in Eq. ( [Til) 
Now, for a polytropic gas. 



(AT) p = Kp'+i, 



(B.2) 



and the equation of hydrostatic equilibrium can be rewrit- 
ten as 



eff- 



This equation integrates to give 
(n + l)p = p($;;^/-$+i|fl xr 



(B.3) 



(B.4) 



where ^^rt is the effective potential at the surface of the 



polytrope. On integrating Eq. ( B.4 ) over the volume of 
the configuration, we obtain 



(n + 1) J p d^r = M$e// -2W+^L-n. 



(B.5) 
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Inserting this last relation in Eq. (p6|), and recalling that 
the pressure vanishes on the surface of a polytrope, we get 



W = 



3A/$ 



(B.6) 



This general expression was derived by Chandrasekhar 
(1961) in a slightly different manner. We now consider 
the case of slowly rotating polytropes. In that case, an ex- 
plicit expression for ^\ff can be deduced from the study of 
Chandrasekhar (1933) on distorted polytropes. The den- 
sity profile can be written 



P = PoO", 



with 



(B.7) 



*• j=i 

The quantities have their usual meaning: po is the central 
density, v = Vi} l^-nGp^^ 9 is the solution of the Lane- 
Emden equation of index n and Pj (/i) are Legendre poly- 
nomials. The functions tpj{^) have been defined and tab- 
ulated by Chandrasekhar (1933). The inner and outer po- 
tentials are given by 



e// 



R 



e 



R 



Co 



(B.IO) 



with R ^ {n + l)KpQ^^^. The constant Aj, Cj are deter- 
mined by requiring the continuity of $ and its derivative 
on a sphere of radius , the first zero of the Emden's func- 
tion with index n. For our purposes, we need to consider 
the equations obtained for j = 0: 



- R^ = Re[ + Rvip'oi^i 



(B.ll) 



(B.12) 



In Eq. (B.ll), we have used ^(^i) = 0. Now, the rela- 
tion between the mass and the central density is given by 
(Chandrasekhar 1933, Eq. (40)): 



M 



-An 



{n + l)K 



AttG 



3/2 



1 



(B.13) 



On the other hand, we define the characteristic radius i?o 
of a rotating polytrope by the relation (see Chandrasekhar 
1933, Eq. (9)): 



Rn 



{n + l)K i_i 



47rG 



1/2 



Inserting Eq. ( |B.12D in Eq. ( |B.13D and using Eq. ( |B.14D 
we obtain 



^Co ^ GM 
£.1 Rn 
On the other hand, we can check that 



Rv = 



(B.15) 



(B.16) 



Substituting these results in Eq. (B.ll), we find that the 
effective potential at the surface of a slowly rotating poly- 
trope is given by 



GM 



2n^Rl 



(B.17) 



We now need to relate the characteristic radius _Ro to 
the total mass M of the system. Eliminating the central 
density po between Eq. ( B.13| ) and Eq. (B.14), and ex- 
pressing V in terms of Rn with the aid of Eq. (B.16), we 
obtain to order f2^: 



Rn — ^0 



1 + if 



SnK- 



271 + 3 2 

2ttG^M-- 



(B.18) 



where Rq is the radius of a non-rotating polytrope. It is 
related to its mass AI by the relation 



(B.9) GNnM'^R^"^ = K. 



(B.19) 



The constant Nn is defined and tabulated in 
Chandrasekhar (1942) for different values of the poly- 
tropic index n. The constant Sn appearing in Eq. (B.18) 
is given by 



1 



3 — 71 V 47r 



1 



(l-™)(3-„) 



- mi 



(B.20) 



and typical values are listed in Table B.l 



From the above results, we can obtain an explicit 
expression for the potential energy of a rotating poly- 
trope. For an axisymmetrical system, we have the relation 
L — I ft where / is the axial moment of inertia (50). To 
our order of approximation, we just need to determine the 
value of / for a non-rotating polytrope. Using the relations 
p = po9" and r — £Ro/£i, expressing the central density 
as a function of Rq by the relation ( |B.14 ) and using the 
mass-radius relation (B.19), we obtain 



/ = InMR, 
with 



•0 7 



47r 



(n-f 1)7V„ 



(B.21) 



(B.22) 



Therefore, according to Eqs. (|B.6|) (|B.17D and (|B.21D , the 
potential energy of a slowly rotating polytrope of index n 
is given by 



(B.14) W 



1 



3GAP 
Rn 



QnMRln^ 



(B.23) 
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Table B.l. The constants of slowly rotating polytropes. 



T 



1 0.63662 

1.5 0.42422 8.73193 

2 0.36475 557.60 

3 0.36394 oo 



0.08320 -0.1009 0.20264 

0.20461 0.40008 0.11127 

0.15485 0.30331 0.05846 

0.07536 0.15087 0.01230 



0.47720 -3.8 lO"'' 0.02257 0.04434 0.00107 



with 

$1 
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where Rq and Rq are given by Eq s. ( |B.19| ) and ( [B.ISD 
respectively. For £1 = 0, Eq. (B.23) reduces to the well- 
known Ritter's formula (see Chandrasekhar 1942). Note 
also that the dimensionless rotation parameter can be 
written 



" GM ■' 



with 



1 

2^ 



47r 
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The polytropic indices n = 1 and n = 3 are special. 
For n = 1, the problem can be solved analytically. The 
following results are well-known (Chandrasekhar 1942, 
1933): 6 = vr, e{ = -I/tt, M^i) = 1, V'(,(6) = 1/vr, 
A^i — 2/7r. According to Eq. (|B.15 ), the radius Rq — 
Rq = (X/27rG')^/^^i is independant of the central den- 
sity and is uniquely determined by the value of K. The 
mass M is arbitrary and independant of the radius. Eqs. 

2 



( |B.2lD ( |B.25| ) remain vahd with Ii = 2{-k'^ - Q)/2,tt'^ and 
Zi — 2/tt'^. The other quantities can be obtained from 
these results (see Table B.l). The case n = 3 corresponds 
to the reverse situation. For $1 = 0, the mass is given by 
Mo = {K/GNsY^^ and the radius is arbitrary. For fl ^ 0, 
we can still consider that the radius i?n = -Ro is a free pa- 
rameter and determine the mass according to Eq. ( B.13D , 
which can be written 



M = M, 



^1/2 n>3 



with 



S" — 
D3 — 



0.03309. 
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We now consider the case of a degenerate Fermi gas at 
zero temperature. As is well-known, this system is equiv- 
alent to a polytrope of index n = 3/2. In addition, the 
constant K is explicitly given by 



5 \47rr/o 



2/3 



(B.29) 



where 770 is the distribution function of a completely de- 
generate Fermi gas. According to the Virial theorem (pi]), 
the total energy E = K + W of this system is E =^-. 
Therefore, using Eq. (B.23), we have 



E 



3GM^ 
' 7Rn 



2o2 



Qy2MR'on 



(B.30) 



According to Eq. (B.l£), the mass-radius relation of a non 
rotating "fermion ball" is 



MRl 



X 



with 



X 



1 



1 



^^3/2 



5.9723 10-3. 



(B.31) 
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125 V47r, 

For a rotating fermion ball we have, according to Eq. 



(B.18) 



Rn — ^0 1 



with 



27rG4M2?72 ; ' 



K= I — I 53/2 ~ 3.9813 IQ-^ 

125V47r/ -^'^ 



(B.33) 
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Then, the energy ( |B.30[ ) can be expressed in terms of M, 
G, rjo and fJ. In order to make the link with the variables 
introduced in Sec. |^, which are normalized by the box 
radius R, we note the relation 



Rq 



(5127r4x)i/V"'/' 



6.6784 



2/3 



(B.35) 



which directly results from Eq. ( |B.3lD and the definition 
of the degeneracy parameter /i. 
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